Local description of the molecular Aharonov-Bohm effect 
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The Aharonov-Bohm effect is one of the most comprehensible examples of quantum non-locality. 
The so called molecular Aharonov-Bohm effect displays great similarities with the latter, but still, 
we show how this effect can be explained using arguments relying solely on locality, whereby we 
mean that the effect can be traced down to a force acting locally on the phase space distribution. 
Our method hinges on studying the system in its momentum representation, and introducing a 
"conjugate gauge potential" which render an everywhere non-zero synthetic magnetic field. The 
resulting Lorenz force induces a transverse current which can be attributed the equivalence of an 
intrinsic spin Hall effect. The idea is demonstrated for the linear E x e Jahn- Teller model and 
applied to the Li3 molecule, for which its corresponding Hamiltonian is obtained by diabatization 
of ab intio determined adiabatic potential energy surfaces. 

PACS numbers: 31.50.Gh, 03.65.Vf 



Introduction. - One of the most well known, and 
also most transparent, examples of quantum non- locality 
is the one by Aharonov and Bohm [l[; An electrically 
charged particle encircling a magnetic solenoid. Despite 
never being in direct contact with the magnetic field, the 
particle acquires a phase shift proportional to the mag- 
netic flux through the loop traversed by the particle. The 
distance between the particle and the solenoid is not of 
importance, solely the fact that the particle's trajectory 
encircles the singular flux line. Non-locality, in general, 
results from the wave nature of the Schrodinger equa- 
tion. In addition to the wave property, this equation is 
manifestly first order in time. Thus, all that is needed 
to achieve the complete system-evolution is the system's 
wave function at some time t — 0, i.e. ip(x.,0). Classi- 
cally, on the other hand, a particle's trajectory depends 
both on its initial position x and its initial momentum p. 
This distinctness of quantum mechanics has its origin in 
position and momentum being conjugate variables, e.g. 
knowing ^(x, t) gives you the corresponding momentum 
wave function <f>(p, t) via a Fourier transform. Conse- 
quently, as the Heisenberg uncertainty relation also tells 
us, affecting one of the two variables will inevitably in- 
fluence the other. 

In this work we demonstrate how a situation seem- 
ingly non-local can be put in a local form by making 
explicit use of the inseparable connection between x and 
p. The system we have in mind is of the Jahn- Teller 
(JT) type @, more precisely the E x e JT model. This 
has been a paradigm model in the chemical and the 
solid state physics communities for more than half a cen- 
tury 0, 3 , and most recently it has also turned out to 
be of relevance in other fields [4J. Within the language 
of molecular gauge theory [j| , an effective non-zero mag- 



netic flux passes through the symmetry point, i.e. the 
conical intersection (CI), of the model. In analog with 
the Aharonov-Bohm effect (ABE) effect, a wave packet 
evolving away from the CI, but still seem affected by the 
flux passing through it has been termed the molecular 
Aharonov-Bohm effect (MABE). Opposite to the ABE, 
Sjoqvist showed how the MABE can be interpreted as a 
local effect by assigning an electric force acting locally on 
the particle- motion Q. We follow a different route and 
introduce a conjugate synthetic magnetic field, such that 
the corresponding Lorentz force is everywhere non-zero 
and give rise to a spin Hall effect |7| . Thereby, put in this 
dual representation, the MABE of linear JT models can 
be readily explained using arguments based on locality. 
It is important to point out that by local we have in mind 
a force acting locally on the phase space distribution. 

Generalized gauge theory of the linear Exe Jahn- Teller 
model. - For the linear Exe JT model, the doubly 
degenerate vibrational mode with frequency lo will couple 
the two degenerate diabatic electronic states according to 

/ P 2 P 2 f) 2 6 2 \ / - 
H EXS = w l^ + ^ + ^ + ^j +k[Q x a x + Q y a y 

where P x , y and Q x . y are two sets of conjugate variables 

with [Q x ,y, Px.y] = i-i an< i k is the JT coupling strength. 
The cr-operators are the standard Pauli matrices; a x = 
|1)(2| + |2)(1|, a y =i(\l){2\ - |2)(1|), a z = |2)(2| - |1)(1|, 
with |1) and |2) denoting the internal diabatic electronic 
states. Throughout this work we use atomic units. 

In standard molecular gauge theory [f|, the x- and y- 
components of the synthetic gauge potential are 



i{^k\d Q M, A y kl) =-i{i> k \dQ v m (2) 
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where the matrix-indices k and I are either 
— or +, and the adiabatic states | ■(/>+) = 
[exp(-i(p/2)\l) -exp(icp/2)\2)]/V2 and = 
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[exp(-i<p/2)\l)+exp(i<p/2)\2)]/V2, with the Q- 
dependent phase ip — arctan (Q y /Q x ). The adiabatic 
states diagonalize the JT coupling part in Eq. ([T]). The 
overall phase ambiguity of the adiabatic states \ip±) 
represents the gauge freedom of the potential ([2]). The 
gauge independent magnetic field can be expressed in 



terms of the field tensor as 
the Levi-Cevita symbol and 
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i[Aj,A k ] 



(3) 



The last part of the field tensor derives from the matrix 
structure of the gauge potential, and whenever the two 
components A x and A y do not commute the system is 
said to be non-abelian. In terms of the E x e JT model, 
this synthetic magnetic field vanishes everywhere except 
at the CI at (Q x ,Q y ) — (0,0), where it becomes singu- 
lar; ergo also the Lorentz force must strictly disappear 
in the punctured Q x Q y -pl&ne K\(0,0). Moreover, since 
the adiabatic states \ip±) are independent on u> and k, it 
follows that also the magnetic field is independent on the 
system parameters. 

Setting the JT coupling strength to zero, k — 0, results 
in a Hamiltonian of identical sets of two-dimensional har- 
monic oscillators. Let x, y, p x , and p y be the quantum 
expectation values of Q x , Q y , P x , and P y respectively. 
For a Gaussian initial state displaced from the origin, the 
trajectories in phase space (x(t),p x (t)) and (y(t),p y (t)) 
will agree with the classical circular ones. To understand 
how the dynamics is affected by the JT coupling, let us 
rewrite the Hamiltonian as 
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where A Ij9 = ——a Xty are the two components of the 
"conjugate" gauge vector potential, and $ = k 2 /cu the 
"conjugate" scalar potential. In this representation, the 
gauge potentials are constant but, due to their non- 
abelian character, give rise to a non-trivial conjugated 
magnetic field and a conjugate Lorentz force, 



k 2 

B z = —i[A x , A y ] = 2—a z 



■^a z {Q y , Q x ). 



(5) 

The interpretation of Eq. Q is a particle with two inter- 
nal states moving in a harmonic potential influenced by 
a constant state-dependent magnetic field; if the particle 
is in the |l)-state it sees a field pointing in the negative 
z-direction and the opposite holds for state |2). Unlike 
the synthetic magnetic field in the position representa- 
tion, the conjugate magnetic field does depend on both 
uj and k. Consequently, it can not be seen as some sort 
of Fourier transform of the latter. Furthermore, put in 
a momentum representation, the JT coupling is nothing 
but a Rashba coupling 8] frequently appearing in con- 
densed matter theories, such as the intrinsic anomalous 
and spin Hall effects Q. 
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Figure 1: (Color online) In (a) the time evolution of the quan- 
tum expectation values (x(t),y(t)) are displayed, while (b) 
shows the corresponding semiclassical averaged trajectories 
obtained using the TWA approach. The final quantum den- 
sity distributions in real and momentum space are depicted 
with contour- lines in (c) and (d), while the dots show the 
semi-classical distributions. 



Numerical results of the linear Exe JT model. - If we 
release a localized wave packet with initial expectation 
values (p x (0),p y (0)) = (0,0) and (x(Q),y{0)) = (ar ,0), it 
appear reasonable to expect y(t) — since the effective 
Lorentz force is strictly zero outside the CI and should 
not affect the dynamics as long as the wave packet is far 
from the CI. On the other hand, the conjugate effective 
Lorentz force of Eq. ([5]) is everywhere non-zero (in mo- 
mentum space) , and the propagation of the wave packet 
must be influenced accordingly. Since the Lorentz force is 
intrinsically acting in the transverse direction of motion, 
the wave packet has to build up a non-zero momentum 
p y {t) in the y-direction. This, in return, will cause a de- 
pletion from the y(t) — trajectory, and our naive guess 
that y{t) = must be wrong. 

We numerically demonstrate this presumption by wave 
packet propagating, using the split-operator method @ 
with an initial Gaussian state, 



i>(Q x ,Q y ,o) 



1/4 



exp 



(Q x - x ) 2 + Q 



2<i 2 



(6) 

with xq = 10 and the width a = 1 corresponding to a co- 
herent state given the frequency of the harmonic oscilla- 
tors. The wave packet is propagated until t = 15 000 a.u., 
and the parameter values, k = 0.01 a.u. and uj = 0.02 
a.u. The resulting trajectories of the expectation val- 
ues (x(t), y{t)) are shown in Fig.[T](a). Directly after the 
propagation has started, a positive y(t) component builds 
up. As soon as the wave packet has crossed into nega- 
tive Q x , the effective force is reversed and y(t) decreases 
again and at the first turning point, x(t) « —10 a.u., it 
has become negative y{t) < 0. When the evolution con- 
tinues, a "bending" of the harmonic motion takes place in 
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which vibrational excitations in the £-mode is swapped 
into excitation in the initially empty y-mode. This can be 
seen as a non-zero transverse current in the ^-direction, 
mimicking the intrinsic spin Hall effect QjJ. The density 
\ip(Qxi Qy, t)\ 2 and its corresponding momentum density 
\cj)(Px , P y , t)\ 2 of the final wave packet are displayed with 
solid contour lines in Fig. [T] (c) and (d) . During the prop- 
agation the wave packet stays localized both in real and 
in momentum space, while for even longer times the wave 
packet, due to anharmonicities originating from the JT 
coupling, will spread out and cover more or less the entire 
accessible phase space 

The transverse current is a result from the JT cou- 
pling having the same structure as a Rashba coupling, 
known to give rise to Hall effects. In order to make the 
phenomenon more transparent we derive a set of semi- 
classical equations of motion, similar to those of Wong de- 
scribing a classical color-charged particle in a non-abelian 
field Qj, 



(7) 



where we recognize s a = (a a ) with a = x, y, z, and the 
dot indicates time derivative. Note that we denote the 
above semi-classical variables the same as the quantum 
expectation values, even though they are not equivalent 
in a strict sense. Indeed, the above set of equations does 
not take higher order quantum correlations into account, 
i.e. the equations are truncated and do not consider 
products like (xs z ), {ys x ) and so on. To extract some 
information about the conjugate Lorentz force, we take 
the derivative of p x (or p y ), p x = —u> 2 p x — k 2 ys z . This is 
the conjugate analog of Newton's second law. The first 
term on the right is the "acceleration" deriving from the 
harmonic motion, while the second is the state depen- 
dent part related to the conjugate Lorentz force being 
proportional to y. 

The semi-classical approach is compared with the full 
quantum wave packet propagation using the truncated 
Wigner approximation (TWA), which has turned out 
to often reproduce the quantum dynamics surprisingly 
well (l3j . Thus, we randomly sample initial values of po- 
sitions and momentum according to the probability dis- 
tributions \ip(Q x ,Q y ,0)\ 2 and \4>(P X , P y ,Q)\ 2 and propa- 
gate them using Eqs. ([7]). The expectation values dis- 
played in Fig. [1] (b) are approximated with the averages 
over all trajectories. The final semi-classical distribu- 
tions obtained from N = 50 000 sampled trajectories are 
depicted in (c) and (d). It is found that this method 
somewhat overestimates the transverse motion, but still 
captures the quantitative features. This discrepancy is 
most likely due to large quantum fluctuations in the spin 
variables which is not captured in the above semi-classical 



analysis. As initial values for the spin isovector we used 
S = (s x ,s y ,s z ) — (0,0,1) for all sampled trajectories, 
i.e. we have not taken initial quantum fluctuations into 
account for the internal two- level system [l4| • 

In addition, for both the quantum and semi-classical 
approaches, we have numerically verified that the trans- 
verse force, i.e. the build-up of a transverse current, 
scales as ~ k 2 for short times as predicted according 
to Eq. flSJ. For longer times, due to the interplay be- 
tween dynamics within the internal states, and dynamics 
arising from the effective Lorenz force and the harmonic 
oscillators, the evolution becomes rather complex, which 
is also known already for free Rashba coupled particles 



for which the problem is analytically solvable 15 1 



It is important to note that for the parameters k and 
u>, the ratio u>/k is larger than that obtained from ab 
initio calculations of the Li3 molecule. As illustrated be- 
low when the dynamics on the two lower states of Li3 is 
explored, the anharmonicity induced by the CI greatly 
destroys the harmonic motion and therefore also the reg- 
ular structure of Fig. [TJ 
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Figure 2: (Color online) The upper plot (a) shows the lower 
energy potential surface of the Li3-trimer (contour lines) to- 
gether with the quantum expectation values x(t) and y(t) 
(thick black line). The lower figures display the distribution 
P(Q y ,t) at time instants, t ~ 410 a.u. (b) and t as 820 a.u. 
(c). The initial state is a Gaussian with y(0) — p x (0) = 
p y (0) — a.u. and x(0) = —0.7 a.u., and width a = 1/2 a.u. 
The final propagation time is approximately 4100 a.u. It is 
particularly clear how the system, as an appearance of the 
spin Hall effect, directly builds up a non-zero displacement in 
the Qy-direction. 
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Effect in Li^-trimer. - In the idealized case of the lin- 
earized E x e JT model, the presence of a transverse force 
is clear from the momentum representation of the prob- 
lem. Once higher order terms are taken into account, 
i.e. the coupling and the potential are no longer lin- 
ear and quadratic respectively, introducing a conjugate 
gauge potential is non-trivial. Nevertheless, we may ex- 
pect that the dynamics should still be greatly influenced 
by effects deriving from linear order terms. Thus, as for 
the above example, the augmentation of a transverse cur- 
rent is taken as a smoking gun for a spin Hall effect which 
in return can be ascribed some synthetic Lorentz force. 

As a demonstrating example we consider the Li3 
molecule which possesses a CI between its two lowest elec- 
tronic states [la]. More precisely, the Li3 molecule pos- 
sesses three vibrational normal modes corresponding to 
bending, asymmetric stretching, and symmetric stretch- 
ing modes, which are all taken into account in our anal- 
ysis. The adiabatic potential energy surfaces are com- 
puted with the MOLPRO program package [TtJ using 
the Multi-Reference Configuration Interaction (MRCI) 
method with molecular orbitals obtained from the state- 
averaged Multi-Configuration Self Consistent Field (MC- 
SCF) calculations with an active space consisting of three 
electrons distributed among fifteen orbitals (all molecu- 
lar orbitals composed of 2s and 2p atomic orbitals). The 
MRCI calculations are performed using the same refer- 
ence space as for the MCSCF calculations. Excitations 
out of the core orbitals as well as single and double exter- 
nal excitations are all included. The calculations are car- 
ried out using the aug-cc-PVTC basis set [l8|- By least 
square fitting of the computed two lowest adiabatic po- 
tential energy surfaces to the eigenvalues of the potential 
part of the JT Hamiltonian, the corresponding JT param- 
eters can be obtained [lj|. Here, in addition to the linear 
terms, the fitted Hamiltonian also contains quadratic and 
cubic anharmonic terms. The obtained JT parameters 
are used to determine the adiabatic to diabatic transfor- 
mation matrix. Then the ab initio computed adiabatic 
potential energy surfaces are transformed into the dia- 
batic potentials and couplings used for the wave packet 
propagation. A reliable description of the system Hamil- 
tonian is then obtained also far from the conical inter- 
section, where the Jahn- Teller Hamiltonian is no longer 
a good approximation (l9| . The details of the numerical 
procedures for the Li3-calculations, together with cor- 
responding results of the potentials, will be published 
elsewhere [201 ] . 

The full three-dimensional wave packet propagation 
is carried out using the Multi-Configurational Time- 
Dependent Hartree method [2l[. The initial state is a 
real valued Gaussian with width a — 1/2 located at 
the lower adiabatic potential surface at (Q x ,Q y ,Qs) — 
(—0.7,0,3.2). For this choice, the initial state is, to a 



good approximation, not in contact with the synthetic 
flux through the CI. The value Q s — 3.2 corresponds to 
a global minima in this coordinate, and the dynamics 
should predominantly appear in the Q x and Q y coordi- 
nates corresponding to bending and asymmetric stretch- 
ing. The lower potential surface is symmetric with re- 
spect to Q y = (see Fig. [5]), and, as for the linear E x s 
JT model, this implies that a non-zero y(t) signals a spin 
Hall effect. The numerical results are presented in Fig. [2] 
showing the trajectory (x(t),y(t)) = {(Q x ), (Q y )) with a 
thick line in (a), and two snap-shots of the projected dis- 
tribution P{Q y , t) = J dQ x \ip(Q x ,Qy,t)\ 2 in (b) and (c). 
The final propagation time in (a) is 4100 a.u. The instan- 
taneous build-up of a transverse force/current is clearly 
demonstrated in the figure, and note that the sign of this 
force is in accordance with our previously introduced con- 
jugated Lorenz force. This is also demonstrated in the 
asymmetry of the projected distributions already seen at 
very short times. Due to large anharmonicities, the wave 
packet rapidly spreads and the effect is therefore only 
visible for short times. 

Conclusion. - Using knowledge of the dual model of 
the Exe JT Hamiltonian, we demonstrated how the cor- 
responding MABE can be attributed only to local forces. 
In the conjugate picture, the inherent gauge potentials 
are non-abelian and constitute an everywhere non-zero 
Lorentz force which causes a wave packet far from the CI 
to deflect from its expected classical trajectory. In con- 
densed matter theories, the same mechanism is essential 
for; the spin Hall effect, topological insulators, high-T 
superconductors, Dirac and majorana fermions, just to 
mention a few examples. In this respect, the present work 
bridge two of the main fields of physics, molecular and 
condensed matter, and our results suggest that known 
effects in one field can be translated into the other. 

Our dual picture analysis is exact in the linear case, 
but is not apparent when non-harmonicity and the extra 
degree-of- freedom are taken into account. Using a full 
quantum calculation, taking all degrees-of-freedom into 
account as well as all higher orders in the JT model, we 
also presented evidences how the intrinsic asymmetry of 
the L13 molecule dynamics can be assigned a spin Hall 
effect. 
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